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Abstract We investigate large-angle scale temperature anisotropy in the Cosmic Microwave 
Background (CMB) with the Wilkinson Microwave Anisotropy Probe (WMAP) data and 
. model the large-angle anomalies as the effect of the CMB quadrupole anisotropies caused 

(-H ' by the local density inhomogeneities. The quadrupole caused by the local density inhomo- 

Q^, geneities is different from the special relativity kinematic quadrupole. If the observer inhabits 

a strong inhomogeneous region, the local quadrupole should not be neglected. We calculate 
such local quadrupole under the assumption that there is a huge density fluctuation field in 
direction (284°, 74°), where the density fluctuation is lO^-^, and its center is 112/i~^ Mpc 
^ I away from us. After removing such mock signals from WMAP data, the power in quadrupole, 

C2, increases from the range (200 ~ 260^K^) to ~ lOOO^K^. The quantity S, which is 
used to estimate the alignment between the quadrupole and the octopole, decreases from 
^ , (0.7 - 0.74) to (0.31 - 0.37), while the model predict that C2 = 1071. S^tK^, S = 0.412. 

• So our local density inhomogeneity model can, in part, explain the WMAP low-^ anomalies. 

I ! 

Q\ ■ Key words: cosmology: cosmic microwave background — cosmology: large-scale structure 

of universe 

(N ' 

o 



1 INTRODUCTION 



Although the WMAP data are regarded as a dramatic confirmation of standard inflationary cosmology 
(Vale 120031 de Oliveira-Costa & Tegmark 2006; Gaztanaga et al. 120031 Gordon et al. 12005), some anoma- 
k> , lous features have emerged (Inoue & Silk , 2006. Campanelli et al. 120061 Dominik et al. .2004 ). Firstly, the 

;_j ■ amplitude of the quadrupole is substantially less than the expectation from the best-fit ACDM standard 

\ model (Abramo et al. l2056l de Oliveira-Costa et al. l2004t Efstathiou 120041), which was found by COBE a 

decade ago (Bennett et al. |1996t and confirmed by WMAP (Spergel et al. i2003i) . Secondly, the quadrupole 
and octopole indicate an unexpectedly high degree of alignment (Spergel et al. 120031 de Oliveira-Costa et 
al. 120041 12OO6I Schwarz |2004l Land et al. |2005l Hansen et al. |2004k : Eriksen et al. |2004l i. 

Recently, many efforts have been devoted to explain the origin of the anomalies. They can be systematic 
error, statistical flukes, improper subtraction of known foreground, or an unexpected foreground (Copi et 
al. 120041 120051 12OO6I 1. The WMAP team claims that there are no unexpected systematic errors (Bennett et 
al. |2003l Finkbeiner |2004| ), and Copi et al. ( |2004ll2005ll2006l ) noted that the anomalies are most unlikely to 
be due to residual foreground contamination. Several authors attempted to explain the anomalies in terms 
of a new foreground (Abramo et al. 120061 Gordon et al. 120051 Bennett et al. l2003l Finkbeiner i2004i Prunet 
et al. 120051 Rakic et al.|2006l l. 

Abramo et al. ( 120061 ) showed circumstantial evidences that an extended foreground near the dipole axis 
could distort the CMB. They proposed that the possible physical mechanism, which can produce such a 
foreground, is the thermal Sunyaev-Zeldovich(SZ) effect. But the SZ model, as presented by them, cannot 
account for the anomalous quadrupole and octopole successfully. Therefore, they thought that the Ress- 
Sciama (RS) effect (Rakic et al. 12006b . or the combination of SZ effect and RS effect may be responsible 
for the foreground. Many other authors suggested that the large-angle anomalies are affected by local inho- 
mogeneities (Tomita |2005b , 12005b ; Vale |2005l l. However, when they applied a model in which the Local 
Group is falling into the center of the Shapley supercluster, the discrepancy between the observed data and 
the model prediction became even worse. (Rakic et al. 120061 Inoue et al. l2007l i. 
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Inoue et al. (2007) explored the large angular scale temperature anisotropies due to homogeneous local 
dust-filled voids in a flat Friedmann-Robertson-Walker universe. They found that a pair of voids with radius 
(2 ~ 3) X 10^h~^Mpc and density contrast 5„i ^ 0.3 might help explain the observed large-angle CMB 
anomalies. While Wu & Fang ( 119941 ) explored the possibility that the CMB is affected by local density 
inhomogeneities basing on Tolman-Bondi model. They calculated the quadrupole amplitude of the local 
collapse model with the general relativity (GR). The results show that the CMB anisotropies from the 
local quadrupole contribution can be different from the special relativity (SR) kinematic quadrupole by a 
factor as large as 3, which depends on the size and density fluctuation of the region the observer inhabits. 
Therefore, if we live in a large density fluctuation area, the local quadrupole might be significant in the 
CMB observations. 

The goal of this paper is to examine whether such local quadruple could account for the observed large- 
angle CMB anomalies in WMAP data. Our analysis is based on the 1-year, 3-year and 5-year WMAP 
Internal Linear Combination maps (Spergel et al. 120061 Hinshaw et al. 12006b (henceforth ILCl, ILC3 and 
ILC5). We try to remove the mock CMB foreground caused by the effect described in Wu & Fang ( 119941 ) 
for each observed CMB map under the assumption that we are in a huge density fluctuation area. The 
parameters of the area we adopted are based on Kocevski & Ebeling (I2006I I and Watkins et al. ( I2008l )'s 
work. We reanalyze the WMAP data by using the multipole vector framework in Section 2. In Section 3, 
we review the estimate of the foreground of Wu & Fang ( 119941 ) and present the result of our examination. 
We conclude in Section 4. 



2 LARGE-ANGLE ANOMALIES OF CMB 

In this section, we re-investigate the anomalies reported from the WMAP maps on very large angular scale. 
As we already remarked, the angular power in quadrupole, C2, is less than expected. To measure C2, we 
expand the temperature anisotropy in terms of spherical harmonics (Campanelli et al. 120061 Copi et al. 120041) 

00 I 

AT(0,,^) = ^ ^ aft„y,™(0,0). (1) 

1=1 m=-i 

And the angular power spectrum is defined as 

m——i 

A simple way to quantify the peculiar alignment of the quadrupole and octopole is to use the multipole 
vectors. In the multipole vector representation, the £ — th multipole of the CMB, Tg, can be written in terms 
of a scalar A(^) and £ unit vectors {{)(^'*) !,■■■,£} (Dominik et al. |2004l Copi et al. |2005l |2006b 

T,«^Wn(t;(^'^).e). (3) 

For the statistical comparison, we use the area vectors 

=^{«) x£,(^>j). (4) 



The alignments between the quadrupole area vector and the three octopole area vectors can be evaluated 
by the magnitudes of the dot products between w'-^'^-^'> and each 

Ai = |to(2;l,2).^(3;l,2)| 

A2 = |t«;(2;1.2).^(3;2,3)| (5) 

The widely used estimator that checks for alignments of the quadrupole and octopole planes is the 
average of the dot products (Dominik et al. 120041 Abramo et al. 120061 Katz & Weeks 12004] Schwarz l2004b 

1 ^ 

S^-J2U (6) 

,:=l 

Given a CMB map, the harmonic components can be easily extracted with the HEALPi)(Q (Gorski et 
al. I2005I I software, and the multipole vectors can be calculated by the code provided by Copi et al. ( |2004| l. 
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Our analysis is based on the 1-year, 3-year and 5-year WMAP full sky maps (ILCl, ILC3, ILC5). The 
values of C2 and S for ILCl, ILC3 and ILC5 are listed in Table 1. We can see that C2 Ues in the range 
(200 ~ 260AiK2). 

In order to compare with the ACDM standard model, 10^ mimic CMB maps are generated with Monte 
Carlos (MC) simulation based on theoretical CMB power spectrum predicted by ACDM model, which is 
gener ated by CAMbO (Lewis et al. I2OOOI 1 pack age with the best-fitting cosmological parameters estimated 
from WMAP (Hinshaw et al. '2009|. In the ACDM model, the power in quadrupole is C2 = 1071.5 ^K^, 
while the power in quadrupole for the ILCl, ILC3 and ILC5 are C2 = 204.4 ^K^, C2 = 260.3 ^iK^, 
C2 = 254.1 /LtK^. Clearly, the WMAP data have a low power in quadrupole compared to ACDM model. 

Fig 1 is the histogram of the S statistics generated from 10^ Gaussian random, statistically isotropic MC 
mock maps. The average value of S from 10^ MC simulations is S'acdm = 0.412, which is much lower 
than the S statistics from WMAP data, that is S* = 0.744 for ILCl, S = 0.700 for ILC3, and S = 0.726 
for ILC5. The final rank in table 1 lists the odds P(5acdm > S) of finding a value among the 10® MC 
maps larger than the one observed, from which one can see that the probabilities are 0.8% for ILCl, 2.1% 
for ILC3, and 1.2% for ILC5. This means that the alignment between quadrupole and octopole for each 
WMAP map is significant. 

These alignments could be explained by an unexpected foreground caused by a local collapse due to 
the second-order effect of the density fluctuation area (Wu & Fang |1994] l. In next section, we will briefly 
discuss this foreground. 



3 HYPOTHETICAL FOREGROUND INDUCED BY SUPER LARGE STRUCTURE 

The CMB temperature anisotropy produced by a locally spherical collapse can be modeled basing on a 
Tolman-Bondi universe solution (Wu & Fang |1994] i. 

Because we are interested in the effect of a local density fluctuation, in the following we only consider 
the case of Xq < X^,, where Xq = x^/te, ccq is the distance between the observer and the center of the 
perturbation, Xc = Xc/te, Xc is the size of the perturbed region. When the initial density perturbation Sq 
is assumed to be constant in the region x < Xc, the first-order solution consists mainly of two parts: a 
monopole term and a dipole term which we are familiar with. The second-order solution of AT/T is (Wu 
& FangUm 
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where To = (1 + Zdf/^ and Zd is the redshift at decoupling time t^, and 4* is the incidence angle of the 
photon. 

When the terms of the order of Sq and Tj|^'^ are taken into account, the quadrupole anisotropy caused 
by local density fluctuation should be (Wu & Fang |1994l l 
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The first term in the left-hand side of equation (8) is the SR kinematic quadrupole anisotropy. Equation 
(8) tells us that if higher orders are involved, the SR kinematic quadrupole may not always be a good 
approximation of the quadrupole produced by a local collapse. The local quadrupole anisotropy strongly 
dependents on the size, matter density in the peculiar field, and the position of the observer Fig 2 shows the 
quadrupole amplitude as a function of the distance between the observer and local gravitational field xq. The 
SR kinematic quadrupole is denoted by solid curve, and the local quadrupole is denoted by dotted curve. 
We assume Xc = 1000/i~^ Mpc to satisfy Xc > xq. The quadrupole showing in Fig. 2 is along the center 
of the perturbation. Fig 3 shows the relationship between the amplitude of local quadrupole and the radius 
of the local gravitational field Xc for xq ~ 112/i~^ Mpc. Xc changes from 15Qh~^ Mpc to 1000/i~^ Mpc. 
Because the distance of the observer to the center of the collapse should at least be greater than the distance 
to the Great Attractor, which is estimated to be 80/igo^ Mpc. Therefore, it would be reasonable to take the 
lower value of Xc — 15Qh^^ Mpc which is about 2 times of the distance to the Great Attractor and the 
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higher value of Xc = 1000/i~^ Mpc which is about the size of horizon (Wu & Fang |1994b . We find that the 
influence of xq on the amplitude of local quadrupole is about one magnitude larger than the influence of Xc- 
When xq is fixed, the results change little with Xc- Fig 4 shows the corrected C2 of ILC5 as a function of 
Xc when xq = I12h~^ Mpc. It turns out that C2 = 1022. for all values of Xc- 

In order to explain the large-angle anomalies we propose a model that we are in a large density fluctu- 
ation area. As Kocevski & Ebeling (I2006I I suggests that 56% of the Local Group's (LG) peculiar velocity 
is induced by more distant overdensities between 130 and 180 Mpc away. Watkins et al. (12008b also notes 
that the bulk flow within a Gaussian window of radius 50 Mpc is 407 ± 81 km s^^ toward / = 287° ± 9°, 
= 8° ± 6°, and roughly 50% of the LG's motion is due to sources at greater depths. Interestingly, we find 
that a region with a density fluctuation S 10^^ over a distance ~ 112/i~^Mpc away on the direction 
of (284°, 74°) may be responsible for the origin of the anomalies on large angular scales. We compute the 
mock foreground (equation (8)) using these parameters. Fig 5 shows the map of the contribution of CMB 
anisotropics caused by the local density fluctuation. 

After subtracting such a mock foreground from the CMB sky maps of the WMAP observation, we find 
that the power in quadrupole will dramatically increase and the alignment of the quadrupole and octopole 
plane will be weakened. In Table 2 we compare the quadrupole and S obtained from the "foreground- 
corrected" WMAP data to those obtained from fiducial ACDM model. The powers in quadrupole of the 
three WMAP maps increase to C2 = 1064.2 ^K^, C2 = 1034.2 /iK^, C2 = 1022.3 ^iK^, respectively, 
which is apparently in much better agreement with the ACDM model. Furthermore, from the S statistics, 
one can see that the frequencies P{Sacdm > S) of finding a ACDM simulation with a S value larger than 
that from WMAP seem to converge to 75.1% for ILCl, 61.5% for ILC3, 62.2% for ILC5. Therefore, if 
such a large scale structure exist, the foreground model presented here can not be neglected. 

We evaluate the probability that the primary quadrupole is cancelled by the local quadrupole. We gen- 
erate 2000 CMB maps, which have random quadrupole orientations, with the HEALPix software, and the 
input theoretical power spectra, C{£), are generated by the CAMB package. Then we combine the fore- 
ground with the random, statistically isotropic CMB maps. We find that about ~ 28% of the quadrupole 
are consistent with the observed WMAP five year values, that is C2 = 223.479 ± 978.36^1 Therefore, our 
model can explain part of the anomalies. But the large errorbar in the quadrupole measurement may also be 
responsible for the large number 28%. 

4 CONCLUSIONS 

In this paper, we have re-investigated the anomalies in WMAP data. The power in quadrupole is found 
to be C2 = 204.4/iK2 for ILCl, C2 = 260.3 ^K^ for ILC3 and C2 254.1 fiK^ for ILC5, while the 
power in quadrupole for the standard ACDM model is C2 = 1071.5 yuK^. It is obvious that the power in 
quadrupole is less than the expected. By comparing the distribution of the S statistics from WMAP data to 
those from 10^ MC simulation mimic CMB maps, we found that they are consistent at the level of 0.8% 
for ILCl, 2.1% for ILC3 and 1.2% for ILC5. These results indicate that the quadrupole and octopole planes 
are aligned strongly. 

We provide a possible explanation for the anomalies in WMAP data by using the foreground model 
caused by a large density fluctuation. The model depends on the matter distribution, and the position of 
the observer. So we assumed that there is a large-scale structure in direction (284°, 74°), the center is 
~ 112h~^ Mpc away from us, and the density fluctuation is 10^"^. After subtracting the mock foreground 
caused by such area from the WMAP data ILCl, ILC3 and ILC5, we found that the power in quadrupole, 
C2, increases to {^^ 1000 ^K^) level, and the S decreases to 0.31 ^ 0.37 level, which agrees with the pre- 
diction from the standard ACDM model. To conclude, the local gravitational collapse might be responsible 
for explaining the origin of the large-angle CMB anisotropy. 

Recently, it has been suggested by many researchers that the local inhomogeneities can account for the 
large angular scales anomalies (Tomita l2005b . 12005b : Vale 12005b . However, none of the proposed models 
can successfully explain the anomalies (Inoue & Silk 12006b . Because it is well known from the GR that 
in a linear approximation, the behavior of a comoving object in an expansion or collapsing metric can 
not be equivalently described as SR Doppler motion if the higher orders are involved. The amplitude of 
the kinematic quadrupole is about 13% of the cosmic quadrupole (Wu & Fang |1994b . Therefore the CMB 
quadrupole anisotropy calculated as an effect of a local density inhomogeneity can not be approximated by 
a SR effect, which is the main reason why we have derived different results from others. 

However, many other specific features of the anomalies have been discovered, such as anomalously 
cold spots on angular scales ~ 10° (Vielva et al. 120041 Cruz et al. 120051) . and asymmetry in the large-angle 
power between opposite hemispheres(Eriksen et al. |2004t Hansen et al., [2004b; Sakai & Inoue 120081 ). We 
have not interpreted these anomalies with our model explicitly, so further research is expected. 
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Table 1 Power in quadrupole C2 and alignments of the CMB maps for ILCl, ILC3 and ILC5. 
The final row shows the expected power in quadrupole and the average value of S statistics of 
the 10^ Gaussian random statistically isotropic CMB maps. P(S) is the probabiUty that a random 
map has a quadrupole-octopole alignment as high as S. 
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Fig. 1 Histogram of S statistics for 10® Gaussian random, statistically isotropic Monte Carlo 
maps. 




Fig. 2 The quadrupole amplitudes as a function of distance between the observer and local 
gravitational field. The solid line indicates the SR kinematic quadrupole, and the dotted line 
represents the local quadrupole. We assume a higher value for X^, that is Xc = lOOO/i^^ Mpc to 
satisfy Xq < Xc- 

Table 2 Power in quadrupole C2 and alignments of the "foreground-corrected" CMB maps for 
ILCl, ILC3 and ILC5. 





C2 (/iK") 


S 


P(S) 


ILCl-corr 


1064.2 


0.317 


75.1% 


ILC3-corr 


1034.2 


0.371 


61.5% 


ILC5-corr 


1022.3 


0.368 


62.2% 


ACDM 


1071.5 


0.412 


50.0% 



Modelling the WMAP low-^ anomalies as an effect of a local density inhomogeneity 



1.6e-05 



1.55e-05 - 



1.5e-05 




300 500 , 700 

[h"^Mpc] 



900 



Fig. 3 The relationship between the amplitudes of local quadrupole and the radius of the local 
gravitational field for xo = 112h^^ Mpc. 
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Fig. 4 Corrected C2 for ILC5, when xq = 112h-^ Mpc. 




Fig. 5 The local quadrupole map. The direction of the local gravitational field is (284° , 74° ), the 
density fluctuation is 10^'^, and it is about 112h^^ Mpc away from us. 



